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Abstract 

Two complementary approaches to the theory of heavy quarkonia 
are discussed. The nonrelativistic potential models give amazingly ac- 
curate predictions, but lack a theoretical justification. The expansion 
in powers of v /c is theoretically very acceptable, but is not as good in 
giving numerical predictions. The importance of combining these two 
approaches is stressed. 



1 INTRODUCTION 

The subject of this presentation are bound states in the bb, bcicb) and cc 
systems. In the old days it was usual to write that the heavy quarkonia 
containing the í-quark, the í-antiquark or both will be the most interesting 
ones . Now, however, it is known that the lifetime of the ŕ-(anti)quark is 
too short for hadronization, so that the quarkonia containing í-(anti)quarks 
do not exist. We limit our discussion to bound states below the thresholds 
for strong decays, i. e. below 2M B = 10.558 GeV for the bb systém, below 
M b + M o = 7.146 GeV for the bc (ôb) systém and below 2M D = 3.690 GeV 
for the cc systém. We also ignore purely relativistic effects like hyperfine 
splittings, fine splittings etc. 
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Plotting the excitation energies M — M( 3 Si) for the three families of 
quarkonia one makes the following two observations. Firstly, the spectrum of 
excitations in the range below the thresholds for strong decays (M th ) depends 
little on the quark masses. Secondly, M t h — Mzg x increases with the quark 
masses. Thus the number of bound states below M th , ignoring the fine and 
hyperfine splittings, is three for čc, four for bc{cb) and ten for bb. Actually, 
many of these states have not yet been observed, but the predictions of the 
nonrelativistic potential models are so reliable that there is little doubt about 
the correctness of this counting. 

We shall discuss two complementary approaches. Nonrelativistic poten- 
tial models give for many quantities predictions in excellent agreement with 
experiment. On the othere hand their relation to accepted theory is unknown 
and they have serious consistency problems. For this reason they are often 
interpreted as recipes rather than as theoretical models. Nevertheless, they 
are generally used, when realistic predicitons are necessary. The recently 
proposed method of expanding in powers of v/c (the heavy quarks are slow) 
[[i]], 0, has much better theoretical foundations. On the other hand its 
predictive power cannot (yet?) compete with that of the nonrelativistic po- 
tential models in the range of their applicability. Nevertheless, this approach 
has produced some generál results of great interest. We shall show how it re- 
solved a problém, which has resisted attempts to clarify it for almost twenty 
years. 



2 NONRELATIVISTIC POTENTIAL MOD- 



There are many nonrelativistic potential models for heavy quarkonia and 
each of them has been introduced for some good reason. Thus e. g., among 
the more than ten published models using potentials of the form 



where the constants a, 6, a, (3 are nonnegative, the Cornell potential a = (3 = 
1 Q becomes Coulombic for r — > and stringy for r —>■ oo as expected from 
QCD; the logarithmic potential of Quigg and Rosner a = [3 —* U gives an 
excitation spectrum, which does not depend on the quark masses; Martin's 
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potential a = 0, (3 = 0.1 makes it possible to perform a particularly 
elegant mathematical analysis. Since the point we want to stress is the pre- 
dictive power of the nonrelativistic potential models, we choose the potential 

i 



V (r) = 0.706380 



0.460442 



+ 8.81715, (2) 



where r -1 and V (r) are in units of GeV. The corresponding mass of the b 
quark is = 4.80303 GeV. This potential and this quark mass have been 
chosen so as to optimize the agreement of the predictions with the following 
measured quantities characterizing the bb systém: the masses M(1S), M(1P), 
M(2S), M(2P), M(3S); the squared moduli of the wave functions at zero 
separation between the quark and the antiquark |^ 1>s (0)| 2 , |^2s(0)| 2 , iV^síO)! 2 ; 
the absolute values of the dipóle electric transition matrix elements | (lP\r\2S) |, 
| (2P\r\3S) | and the ratio of such absolute values | (lS\r \2P) |/| (2S\r\2P) |. Ac- 
tually, instead of fitting the eleven observables with four free parameters, we 
fitted eight observables with one partameter and fixed the other three 
parameters so as to reproduce exactly the remaining three observables. This 
has little effect on the quality of the fit and makes the analysis simpler and 
more meaningful statistically. Another advantage of this procedúre is that 
the poorly known corrections in the evaluation of the squares of the wave 
functions at the origin from the directly measured leptonic decay widths 
almost drop out. 

For the parameters quoted above one finds x 2 = 6.5 for seven degrees of 
freedom, which means an excellent fit. This is nontrivial, because the ac- 
curacy of the experimental data is high. In particular, the four masses are 
measured with uncertainties ranging from 0.2 MeV to 0.5 MeV, i. e. with an 
error below 0.01%!. Incidentally, our parmeters are given with six digit accu- 
racy, only in order to enable the reader to check our computer calculations. 

One of the early successes of the nonrelativistic quark models was that 
the potentials tuned to describe the 5c quarkonia have also given an accept- 
able description of the bb quarkonia, when these were discovered several years 
later. This fact, which had been expected in the nonrelativistic quark models 
"QCD is flavour blind", is not predicted any more in the modern approach 
based on the expansion in powers of velocity. Therefore, it is of interest to 
check, how this universality is satisfied with the present data. Introducing 
one more parameter m c = 1.3959 GeV and assuming that the constant in 



3 



the potential changes by 2(m c — m&), we have calculated for the cc quarkonia 
the values of the observables M(1S), M(LP), M(2S), |^is(0)| 2 , \*p2s(0)\ 2 
| (líSlr 1 1-P) | and |(lP|r|2S')|. The masses agree with experiment within 4 
MeV. This margin is about an order of magnitude more than the experimen- 
tal uncertainties, but on the absolute scale it is not much. The wave functions 
and the electric dipóle matrix elements have to be rescaled (a rescaling fac- 
tor of 1.3 for the wave functions and a rescaling factor of 0.73 for the matrix 
elements), but then they agree with experiment. Some rescaling is expected, 
because when extracting the quantities to be compared with the model from 
the directly measurable ones we used the same correction factors as for the 
bb systém, which in the case of the wave functions is certainly, and in the 
case of the matix elements very probably, a rather crude approximation only. 
We conclude that deviations from the bb — cc universality are clearly seen, 
but that they are not very large. 

Let us conclude by stressing that the nonrelativistic potential model can- 
not be interpreted naľvely as a consistent theory. E. g. within this model 
one can calculate the kinetic energies of the quarks and from that their 
mean square velocities. The results are (v 2 ) « 0.25c 2 for the ground state 
of the cc systém and (v 2 ) ~ 0.08c 2 for the ground state of the bb systém. 
This is inconsistent with the assumption that high precision results can be 
derived assuming that the motion is nonrelativistic. To be sure, one can 
speculate that relativistic corrections can be absorbed into redeŕmitions of 
the nonrelativistic parameters of the model, but before this is demonstrated, 
the main argument in favour of the nonrelativistic potential models remains 
their amazing success in predicting experimental results. 

Let us consider now a more formal approach. 

3 EXPANSION IN POWERS OF VELOC- 
ITY 

Since confmement is a nonperturbative effect, there is little hope of obtain- 
ing a good description of heavy quarkonia by summing Feynman diagrams. 
Another obvious idea - to try an expansion in inverse powers of the heavy 
quark mass - is also unlikely to work. Arguments derived from quantum field 
theory can be found in ref 0. Here we shall qualitatively discuss the high 
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mass limit using the Cornell potential 



V{r) = ar h c, (3) 

r 

where a > 0, b > 0. It is plausible that in the high mass limit the systém 
becomes Coulombic. A Coulombic systém with particle masses M has a 
rádius of order i , so that in the high mass limit it is consistent to neglect the 
ar term. Let us try, however, to consider this term as a perturbation. For a < 
the potential is unbounded from below. There is no ground state and the 
perturbation serieš cannot converge. Theorems about the convergence radii 
of power serieš guarantee that if the perturbation serieš in a is divergent for 
some a < 0, it will also diverge for all positive a from the range < a < —a . 
Thus a = is a singulár point. This could be harmless. E. g. the perturbative 
expansion in QCD is believed to be an (asymptotic) power serieš expansion 
around a singulár point. The reason for its success is that the first terms of 
this expansion approach the correct result so rapidly that they are enough 
for most practical applications and, therefore, the convergence or divergence 
of the whole serieš is irrelevant. In the quarkonium problém, however, the 
high mass limit - with the wave function localized in an infinitesimal región 
around r = - is so remote from any plausible description of a quarkonium 
that there is little hope that the first few terms of the expansion will make 
it realistic. Thus, one has to look for another idea. 

Let us note an important implication J^] of the observation that in the 
high mass limit the quarkonium is approximately a Coulombic systém. Since 
the kinetic and the potential energies should be comparable, we expect for 
high masses 

« Mv\ (4) 
R w 

where R = ^- is an estimate of the quarkonium rádius. Since a s (//) is a 
decreasing function of its argument and since v (which is in units of the 
velocity of light c) is less than one, this implies 

v > a. (M). (5) 

Thus, it is inconsistent to include radiative corrections of order 0(a k s ) without 
also including the relativistic corrections of order at least v k . 
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The new and promising idea , @ , is to combine factorization and an 
expansion in powers of the quark velocity (in the rest frame of the quarkonium 
and in units of c) v. Thus, any probability amplitúde, e. g. the probability 
amplitúde for the decay of quarkonium QQ into a gluon pair, is represented 
as a sum of terms. Each of these terms is a product of a soft matrix element, 
which in principle is obtainable from a lattice calculation, but at present is not 
known, and of a hard matrix element, which can be calculated perturbatively. 
Using suitable scaling rules it is possible to ascribe to each term an order n, 
which means that for v — > this term is of order 0(v n ). At each order there 
is a finite number of terms only. The leading term approximation consists of 
all the terms of lowest order. It may be systematically improved by including 
higher and higher order terms. 

In this approach the soft matrix element depends on the quark mass 
in a way, which is beyond our control. Thus, the simple universality from 
the nonrelativistic quark model is lost. Nevertheless, many approximations 
known from other approaches, like spin symmetry, vacuum saturation or 
some of the relations between the production and decay amplitudes can in 
many cases be rigorously justified at sufficiently low orders of the expansion in 
powers of v. Moreover, the results have often simple physical interpretations. 

Perhaps the most impressive success of this approach is the analysis of 
the decays of P-wave quarkonia into light hadrons [H], [j]. This process had 
been studied for a long time in the framework of the nonrelativistic potential 
models. One finds that, since for P-states (^(O)] 2 = 0, Q and Q have small 
probability to get so close to each other as to be able to annihilate with a sig- 
nificant probability. As a result, the decay probability amplitúde is reduced 
(compared to the decay probabilities of the S'-states) by a factor of order 
O (v). The problém is, however, that this comparatively small amplitúde has 
infmite QCD corrections! The infrared divergences in the calculation do not 
cancel ||. 

In the spirit of the expansion in powers of v, we must look for other 
contributions to the decay amplitúde, which are of the same order in v. 
There is one more such term, where the quarkonium component consists 
of a QQ systém in a colour octet S'-state accompanied by a "dynamical" 
gluon so that the whole quarkonium is a colour singlet as it should. This 
component is small - of order 0(v) as compared to the main term. Its QQ 
part, however, being an S'-state annihilates easily, so that the contribution of 
this component to the decay amplitúde of the quarkonium into light hadrons 
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is of the same order in v as the decay amplitúde of the P-wave term. Thus, 
both terms must be included in a correct leading order calculation. After 
this is done, the infrared singularities in the QCD corrections cancel || and 
a finite, acceptable result is obtained. Thus, at least in principle, the problém 
posed in 1976 has been solved. 

4 CONCLUSIONS 

The expansion in powers of the velocity f is a respectable physical approach. 
It can be used to prove statments, which in the nonrelativistic potential ap- 
proach had to be, often implcitly, conjectured. It also introduces important 
corrections beyond the scope of the potential approach, like the contribution 
of the | (QQ) , gluon) state to the annihilation of the P-wave quarkonia into 
light hadrons. Nevertheless, for making practically useful predictions it can- 
not, at present, replace the potential models. The key problém seems to be, 
how to combine the advantages of the two approaches? How to calculate 
the soft matrix elements of the expansion in powers of the velocity v from 
potential models, or how to justify the potential models using the powerful 
formalism of the expansion in powers of v? 
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